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Search for the proton decay mode p —>• uK + with KamLAND 
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We present a search for the proton decay mode p -A VK + based on an exposure of 8.97 kton-years in the 
KamLAND experiment. The liquid scintillator detector is sensitive to successive signals from p —> VK + with 
unique kinematics, which allow us to achieve a detection efficiency of 44%, higher than previous searches in 
water Cherenkov detectors. We find no evidence of proton decays for this mode. The expected background, 
which is dominated by atmospheric neutrinos, is 0.9 ± 0.2 events. The nonbackground-subtracted limit on the 
partial proton lifetime is r/B(p -A vK + ) > 5.4 x 10 32 years at 90% C.L. 

PACS numbers: 12.10.Dm, 13.30.-a, 12.60.Jv, 11.30.Fs, 29.40.Mc 


I. INTRODUCTION 

Extensions of the Standard Model (SM) of particle physics, 
such as supersymmetric grand unified theories (SUSY- 
GUTs) EL [i], predict that the three SM coupling constants 
are unified at a scale of ~2 x 10 16 GeV. Proton decay is 
an indispensable probe to determine the unification scale ex¬ 
perimentally. In SUSY-GUT models, the proton decay mode 
p -A 17K + is expected to be the leading process through 
color-triplet Higgs exchange in dimension-five operator in¬ 
teractions 0 . At the same time, the predicted decay rate of 
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p -A e + 7 r° can be highly suppressed by the masses of the de¬ 
cay daughters, and is consistent with the current experimental 
lower limits on the partial lifetime M. 

Previous searches for the p -A 17K 3 decay mode were per¬ 
formed by the Super-Kamiokande Collaboration with a fidu¬ 
cial mass of 22.5 kton ffl. They have reported the most 
stringent limit of r/ B(p -A 17K + ) > 5.9 x 10 33 years (90% 
C.L.) with a 260 kton-year exposure 0. In water Cherenkov 
detectors, the direct signal from the K + is not visible because 
its momentum is below the Cherenkov threshold, and the in¬ 
formation on the energy released by the proton decay event 
is lost. So the Super-Kamiokande search relied on three de¬ 
tection channels in which the K + is identified via its decay 
daughters: (1) K + -A in delayed coincidence with 

deexcitation 7 rays from the remnant nucleus from p decay 
within an oxygen nucleus; (2) K + -A li + iy,, monoenergetic 
muon search; and (3) K + -A 7 r + 7 r° search. The search with 
method ( 2 ) is less sensitive than the other methods due to the 
large atmospheric neutrino background. The background-free 
methods (1) and (3) show low efficiencies of 9.1% and 10.0%, 
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respectively, due to backgrounds around the deexcitation 7 ray 
energy and the limited resolution for reconstructed pion mo¬ 
mentum. The systematic uncertainties on the efficiencies are 
22% (Method 1), mainly arising from the dependence on nu¬ 
clear models to predict the probabilities of the deexcitation 7 
rays of the remnant nucleus, and 9.5% (Method 3) from the 
uncertainties of the 7 r + interaction cross section in water and 
the event reconstruction. 

Liquid scintillator (LS) detectors, on the other hand, are 
expected to achieve higher detection efficiency with a lower 
systematic uncertainty because the scintillation light from the 
K + can be readily observed, resulting in a clear delayed coin¬ 
cidence signature with the signals from the K + decay daugh¬ 
ters. In addition, the kinetic energy distribution of the decay 
K + measured by scintillation will peak at ~105 MeV, with a 
spread due to nuclear effects and detector resolution, and will 
provide strong evidence to claim the p —> VK + detection. 
The potential of a large volume LS detector was discussed in 
Refs. Gum, and the detection efficiency of the proton decay 
mode p —> VK + is expected to reach 65% 0 . However, this 
technique has not yet been demonstrated experimentally. 

This paper demonstrates the principle of the new search 
method by presenting the first experimental search results 
for the proton decay mode p —> vK + using a large-volume 
LS detector, KamLAND. KamLAND is the largest scintil¬ 
lation detector constructed to date, with a detector mass of 
1 kton 0 . As discussed in Ref. iflTll , a significant challenge 
is that time difference between the prompt K + and delayed 
decay daughter signals on the same time scale as the scintilla¬ 
tion response, requiring that the two signals be reconstructed 
from within an overlapping detected light profile. We show 
that this challenge can be overcome in practice. The sig¬ 
nal efficiency and background studies presented in this paper 
are also important for a realistic estimate of the search sen¬ 
sitivities in future large LS detectors, and will provide useful 
knowledge for a comparison of the sensitivity between water 
Cherenkov and LS detectors in the future. 


II. KAMLAND DETECTOR 

KamLAND is situated 1 km (2700 m water equivalent) un¬ 
derground in the Kamioka mine in Japan. The detector con¬ 
sists of an 18-m-diameter stainless steel spherical shell which 
defines the boundary between the inner and outer detectors 
(ID and OD respectively). The inner surface is instrumented 
with 1325 fast-timing 17-inch and 554 20-inch Hamamatsu 
PMTs facing inwards to a 13-m-diameter volume of LS. The 
density of the LS, which consists of 80% dodecane, 20% 
pseudocumene (1,2,4-trimethylbenzene), and 1.36 g/liter PPO 
(2,5-diphenyloxazole), is 0.780g/cm 3 at 11.5°C. The LS is 
contained in a 135-//m-thicl< transparent nylon/EVOH (ethy¬ 
lene vinyl alcohol copolymer) balloon, which is mechanically 
supported by Kevlar ropes. Outside the balloon is a buffer 
region of dodecane and isoparaffin oils that shields the target 
from radioactive decays emanating from the walls of the de¬ 
tector and the PMTs. The OD is a cylindrical region of pure 
water instrumented with 225 20-inch PMTs. The OD tags 


cosmic-ray muons entering the detector and also absorbs neu¬ 
trons from the surrounding rock. 

The 1325 17-inch PMTs alone give a photocathode cov¬ 
erage of 22%. The total PMT photocathode coverage of the 
inner detector is 34%. Only the fast timing 17-inch PMTs 
were used in this analysis. When a PMT channel registers a 
“hit,” the charge as a function of time is recorded by an analog 
transient waveform digitizer (ATWD) chip. The KamLAND 
ATWDs each capture a 128 sample waveform (about 200 ns). 
Each channel has two sets of three ATWDs, which record 
waveforms at three available gain levels, providing enough 
dynamic range to record signals between one and thousands 
of photoelectrons (p.e.’s). The dual (2 x 3) ATWD system al¬ 
lows signals to be captured on a second set of ATWDs during 
the 30 ps required to read out the recorded waveforms from 
the first light detected in an event. 

The data reported here are based on a total live time of 
9.82 years, collected between lanuary 10, 2003 and May 1, 
2014. The data set is divided into two periods: period I 
(3.79 years live time) refers to data taken up to May 2007, at 
which time we embarked on a LS purification campaign that 
continued into 2009. Period II (6.02 years live time) refers 
to data taken during and after the LS purification campaign. 
We found that the scintillation and optical parameters are dif¬ 
ferent between the two periods, so we analyzed the data for 
each period separately, and finally the results are combined 
to estimate the partial proton lifetime. We use a LS target 
mass of 913 tons for the p —? vK + search, resulting in a 
total exposure of 8.97 kton-years. Since September 2011, 
the KamLAND-Zen experiment was launched to search for 
neutrinoless double beta decay using a newly introduced (3(3 
source, 13 tons of Xe-loaded LS in a 3.08-m-diameter inner 
balloon at the center of the detector 0 . The volume fraction 
of Xe-LS to all LS is 0.014 and the concentration of enriched 
xenon gas in Xe-LS is less than 3 wt%, so the impact of the 
presence of the (3/3 source on the proton decay analysis is neg¬ 
ligible. 


III. SIMULATION 

Simulations were performed in order to understand the 
event characteristics from p —> vK + signals and backgrounds 
in KamLAND. As described below, we use Monte Carlo 
(MC) simulation tools, an event generator for high-energy 
event studies, and a general particle tracker applied to the 
KamLAND detector. In the p —>- VK + search, the remaining 
background after the cosmic-ray muon rejection is dominated 
by atmospheric neutrinos. The atmospheric neutrino back¬ 
grounds need to be sufficiently removed by requiring delayed 
coincidence cuts based on the scintillation signals from the 
prompt K + and the delayed decay daughters, while maintain¬ 
ing a certain level of the detection efficiency. The optimum 
selection was studied by comparing the event characteristics 
of MC samples from p —>■ VK ! decays and atmospheric neu¬ 
trinos. Finally, the signal efficiency and the background rates 
were evaluated from MC samples. 
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FIG. 1: The K + kinetic energy distribution from bound proton de¬ 
cays in 12 C with (solid line) and without (dashed line) the intranu¬ 
clear cascade. The monochromatic line indicates the K + kinetic en¬ 
ergy for free proton decays. 


A. p —> vK + decay 

The signal efficiencies for p —> vK + decay in KamLAND 
were estimated with MC simulations. We generated 10,000 
event samples distributed in the complete LS volume (1171 ± 
25 m 3 ). For a free proton at rest, this two-body decay results 
in the back-to-back emission of the V and K + with monochro¬ 
matic kinetic energies of 339 MeV and 105 MeV, respectively. 
However, the prompt signal from bound proton decay in 12 C 
will be modified by various nuclear effects, namely the pro¬ 
ton’s Fermi motion, the nuclear binding energy, and intranu¬ 
clear cascades triggered by Ff + -nucleon interactions. In order 
to model these effects, we used a nucleon decay event genera¬ 
tor in the GENIE-based MC code (version 2.8.0) 10. GENIE 
is usually used for neutrino interactions with nuclear targets. 
However, the code provides models of nucleon states within a 
nucleus as well as intranuclear hadron transport after neutrino 
interactions, which are essential also for nucleon decay sim¬ 
ulations. Previous MC studies in Ref. 0 were performed 
under the condition that the emitted K + from proton decay 
has a kinetic energy of 105 MeV without nuclear effects. 

The distribution of momenta and binding energies of the 
nucleons in 12 C is modeled by the spectral function of Ben- 
har 0 , which reproduces the experimental data of electron¬ 
ic scattering. In the original GENIE code, the two-body de¬ 
cay of a nucleon was calculated based on a simple Fermi gas 
model ignoring the shell structure of the nucleus and nucleon- 
nucleon correlations, and the modification of the decay nu¬ 
cleon mass by the nuclear binding energy was absent. We im¬ 
plemented the mass correction m' N = uin — Eb (where Eb is 
nuclear binding energy) and a corresponding Fermi momen¬ 
tum based on the spectral function to provide a more realistic 
description of nucleon states in 12 C. 

The produced hadron from nucleon decay is simulated by 
the intranuclear cascade model in GENIE. It describes the 
hadronic interactions with each nucleon inside the remnant 


nucleus before exiting. A bound proton decay of p —> vK + 
in 12 C produces K ’ with low momentum, below 600MeV/c, 
so only elastic scattering ( K + p —>• K + p or K + n —> K + n ) 
and charge exchange ( K + n —> K°p ) with nucleons in 11 B 
can occur. In order to describe the K + interaction with nu¬ 
cleons more accurately, we modified the GENIE code to im¬ 
plement elastic scattering and charge exchange reactions cal¬ 
culated using the data-driven cross sections from partial-wave 
analysis 111611 171. The K + loss probability due to charge ex¬ 
change is calculated to be only 2.1%. Figure D] shows the re¬ 
sulting K + kinetic energy distribution for bound proton de¬ 
cays in 12 C. The smearing of the peak is mainly due to the 
Fermi momentum. For the case of /\ + -nucleon interactions, 
we calculated the energy loss of K + and the energy transfer 
to secondaries, such as protons or neutrons, which have a dif¬ 
ferent scintillation response. 

The initial K + generation with MC is followed by a set of 
processes for particle tracking, scintillation photon emission, 
and photon propagation using a separate MC code based on 
Geant4 1181 . The detector response MC code was tuned to 
reproduce the data of radioactive source 7 rays for low ener¬ 
gies and atmospheric neutrinos for high energies. The scin¬ 
tillation and optical parameters for the MC, including Birk’s 
constant, the light yield, the attenuation length, and reemis¬ 
sion and scattering probabilities, are determined separately for 
period I and period II. After the MC tuning, we found the 
fractional difference of the average detected number of p.e. 
between the measured data and the MC for central-axis ra¬ 
dioactive source calibrations with 60 Co to be less than 2%, as 
shown in Fig. [2] The time jitter and waveform for each gain 
channel for the 17-inch PMTs are measured with dye-laser 
calibration data. These are used to simulate the photon arrival 
time in the detector response MC output. 



FIG. 2: Fractional difference of the average total detected p.e. of 
60 Co events relative to the detector center for data (squares) and MC 
(circles). The decrease of charge around the ±6 m region is due to 
the absence of PMTs in the detector entrance at the top and bottom 
region, and shadowing by suspension ropes. 
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B. Atmospheric neutrino background 

The backgrounds from atmospheric neutrinos are estimated 
with MC simulations. We used the GENIE-based MC event 
generator assuming the atmospheric neutrino flux of Ref. [19] . 
Three-flavor neutrino oscillations based on the best fit param¬ 
eters given in Ref. Eoll were taken into account. The de¬ 
tector simulation after the event generation is performed by 
the Geant4-based MC code described above. We produced 
20,000 atmospheric neutrino MC events in KamLAND which 
corresponds to about 100 years live time. 

The reactions are dominated by quasielastic scattering be¬ 
low 1 GeV, producing a charged lepton in a charged-current 
interaction or a recoil nucleon in neutral current scattering. 
While those reactions overlap with the p —> vK + signal in 
visible energy, the background can be effectively suppressed 
by the delayed coincidence requirement. The background rate 
from the misidentification of single-pulse atmospheric neu¬ 
trino hit patterns as multipulse events is estimated in the fol¬ 
lowing section. 

Other potential background sources are kaons produced 
by the neutrino interactions. The GENIE code simulates 
kaon production via resonances based on the Rein-Sehgal 
model ®. These reactions conserve strangeness, resulting 
in pair production of kaons and hyperons in the final state, 
such as Vp + n — >• p~ + A + K + . Based on the GENIE cal¬ 
culation, the number of events below 2 GeV with a K + pro¬ 
duced via a resonance is expected to be only 0.2 in the data 
set. On the other hand, nonresonant kaon production, such as 
+ p —» p~ + p + K + , is not considered in the GENIE 
code. Although experimental confirmation of such reactions 
is absent, these cross sections are calculable theoretically [22]. 
While they are comparable with the cross sections for resonant 
K + production, due to the lower reaction thresholds, nonres¬ 



onant kaon production could be a more important background 
in the p —> VK + search. However, based on the cross sections 
in Ref. El, the number of events below 2 GeV is expected to 
be only 0.14 in the data set. Therefore, the background con¬ 
tributions from both resonant and nonresonant kaon produc¬ 
tion are not significant for the current p —t VK + search in 
Kam L AND. 


IV. ANALYSIS 

Clear identification of the two scintillation signals from the 
prompt K + and the delayed decay daughters requires the mul¬ 
tipulse fit to the PMT hit-time pattern as discussed later in this 
section. However, such fits to all the recorded events take a 
lot of computing time. Therefore, we apply three selection 
criteria prior to the delayed coincidence fitting: (i) the recon¬ 
structed energy ( E ) based on the PMT charge recorded within 
a restricted 150 ns event window is required to be more than 
30MeV; (ii) the number of OD PMT hits is less than 5; and 
(iii) the small contribution of events in time coincidence with 
the neutrino beam of the T2K experiment El are rejected. 
Unlike the low-energy radioactive decays that comprise the 
bulk of the KamLAND data, most of the p —» Vl\ + decays are 
expected to deposit a much higher energy, typically peaked at 
260 MeV including short-lived decay daughter contributions. 
So the radioactivity background can be effectively suppressed 
with (i). The event rate of cosmic-ray muons which pass 
through the ID is about 0.3 Hz, but those muons can be re¬ 
jected by the OD tagging of (ii). The inefficiency of the OD is 
estimated to be about 0.2%, and the remaining muons are re¬ 
jected by the PMT hit-time distribution, discussed later. The 
events remaining after application of these selection criteria 
are dominated by atmospheric neutrinos. Figure [3] shows the 
reconstructed energy spectrum of candidate events within a 
5-m-radius volume after the prior selections. The observed 
spectrum is consistent with the expectation from the MC sim¬ 
ulation. The uncertainty of the energy scale model, discussed 
in Sec.[V] is constrained from these data. 

These high-energy candidates were used for the p —t 17 K + 
signal search. In order to identify the successive signals from 
K + and its decay daughters over a short time, we decompose 
one event into multiple pulses based on a binned-^; 2 fit to the 
sum of the waveforms from all hit 17-inch PMTs. The wave¬ 
form is a function of the photon arrival time after subtracting 
the time of flight (TOF) from the vertex to each PMT. The 
average waveform for high-energy candidates is used as the 
reference pulse for the fit. Considering the energy and posi¬ 
tion dependence of the waveform shape, the reference pulse 
was prepared for each energy-position bin. The binned -\ 2 for 
M pulses is defined as 


FIG. 3: Reconstructed energy spectrum of atmospheric neutrino can¬ 
didates (black dots) within a 5-m-radius volume. The atmospheric 
neutrino MC (solid line) is normalized to the live time of data. The 
neutrino interaction uncertainties evaluated by the GENIE code are 
indicated by the shaded region. The expectations for the linear en¬ 
ergy scaling with ±1 <t are also shown (dashed lines). 
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At 2 ^), (1) 


where / (T)) is the observed waveform sum for the 1th bin 
(bin width is 1 ns) as a function of T (0 < T < 150 ns), 
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FIG. 4: Typical waveform shape (shaded) expected from (a) and (b)p —► VK + MC example events and (c) an atmospheric neutrino MC event, 
together with the best fit curves from the multipulse fit (solid thick line) and from the single-pulse fit (solid thin line). In the upper panel, the 
residuals from the best fit are shown to identify multipulse events. The fitted energies and, for the multipulse fits, time differences between the 
first pulse from the K + (red dashed line) and the second pulse from the /i + decay daughter of the K + (blue dashed line) are also shown. The 
shape parameters (see text) for each fit are (a) log (L s h ap e) = 0.06 and (b) log (L s ha P e) = 0.02. 





FIG. 5: Typical waveform shape (shaded) from: (a) and (b) muon + decay electron data and (c) atmospheric neutrino data, together with the 
best fit curves from the multipulse fit (solid thick line) and from the single-pulse fit (solid thin line). In the upper panel, the residuals from the 
best fit are shown to identify multipulse events. The fitted energies and, for the multipulse fits, time differences between the first pulse from 
muons (red dashed line) and the second pulse from decay electrons (blue dashed line) are shown. 


is the reference pulse, E is the total visible energy estimated 
from the observed charge within the 150 ns event window, Ej 
is the energy contribution for the yth pulse, and A Tj is the 
time difference between the first and the yth pulse (ATi = 0). 
ct 2 (Tj) is calculated by 


o*(Tj = 


M r 

E 

3 


El 

E 


A - A Tj) 


( 2 ) 


where A <j>(Ti) is the statistical fluctuation (ler deviation) of 
the reference pulse evaluated by the MC simulation. Since at¬ 
mospheric neutrino interactions can produce electromagnetic 
showers, the waveform shape may change. In order to ac¬ 
count for such fluctuations, the parameter a is introduced to 
scale the time (T —>• aT) of the first pulse and is floated to fit 
the peak time within the uncertainty. In the main branching 


ratio of the K + decay channel, K + —> jJi + v^ (63.55%), the 
event can contain up to three pulses, where the second pulse 
accounts for the energy deposition of the /i + , and the third 
pulse is for the e + decay daughter of the /i + . However, due 
to the short digitization window (a few hundred ns) relative to 
the muon lifetime (2.2 /is), and due to ATWD dead time and 
misfires resulting from baseline instability in very high PMT 
charge events, the separated e + event is difficult to identify. 
So in this analysis the e + tagging is not used and we restrict 
the multipulse fit to M = 1 or 2 only. In the second-largest 
branching ratio, K + —>• 7r + 7r° (20.66%), the second pulse is 
for 7T + 7r° decay daughters of the K + with the instantaneous 
electromagnetic 7r° decay, resulting in a higher energy deposit 
than the main branch. The subsequent /i + decay daughter of 
the 7r + has a negligible impact on the fit because of the small 
energy deposit of 4.2 MeV. Similarly, we use the multipulse 
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TABLE I: Main decay modes of K + in descending order of branch¬ 
ing ratio I20h . The prompt K + signal is followed by delayed signals 
from decay daughters of the K + . Their kinetic energies which con¬ 
tribute to instantaneous scintillation are summed up. 


Decay mode 

Branching ratio (%) 

Kinetic energy sum (MeV) 

K + — >■ p + Vfi 

63.55 ±0.11 

152 

K + —> 7r + 7r 0 

20.66 ± 0.08 

354 

A' + — > n + TT + TV~ 

5.59 ±0.04 

75 

K + —> -K°e + v e 

5.07 ±0.04 

265-493 

K + —> iv 0 p + u ll 

3.353 ±0.034 

200-388 

A' + —> Tv + n°Tv° 

1.761 ±0.022 

354 


fit to search for all the K + decay channels. The strongest de¬ 
cay branches and the corresponding kinetic energy sums that 
comprise the delayed scintillation signal are listed in Table[T| 
Figures 0 a) and 0b) show multipulse fit results for typical 
p —>• VK + MC events. The plots are labeled with the best 
fit energies of the prompt K + and delayed /j + scintillation 
pulses (£ 12 ), and the time difference between them (AT). 
To demonstrate the performance of the fit, Figs. 0a) and 0b) 
show multipulse fit results for typical multipulse events from 
muon decay events recorded by KamLAND, labeled with best 
fit values for the reconstructed muon and Michel electron en¬ 
ergies (Ei and E 2 , respectively) and the decay time (AT). We 
also show the performance of the fit on typical atmospheric 
neutrino events in both MC and data in Figs. 0c) and 0c), 
which are best fit with M = 1. 

We require the following criteria to select the p —» VK + 
signal: 

(1) The distance from the center of the detector to the recon¬ 
structed vertex ( R ) is less than 6.5 m, which corresponds 
to the full LS volume. 

(2) The PMT hit-time distribution is required to be more 
pointlike than tracklike to reject cosmic-ray muons which 
have long tracks in the ID. Specifically, the time-based 
vertex-x 2 (Xt), which compares the observed PMT hit¬ 
time distribution as a function of distance to the event to 
the expected distribution, is required to be less than 1.2 
(pointlike) for the full volume, or to lie between 1.2 and 
5.0 (short-tracklike) for z < 5.0 m and x 2 + y 2 < 25 m 2 . 

(3) The best fit time difference between the first pulse and the 
second pulse (AT) is between 7 and 100 ns, considering 
the K + lifetime of 12.4 ns. 

(4) The best fit energy of the first pulse (£j) is between 30 
and 300 MeV, the second pulse energy (£ 2 ) is between 70 
and 600 MeV, as constrained by the decay kinematics. 

(5) The waveform shape is multipulselike to reject single¬ 
pulse events. We use a waveform shape parame¬ 
ter defined as log (L s h ape ) = log(Ax' 2 ) — A, where 
A X 2 = xl,min - xlmin ( where XM,min is the beSt fit X 2 



0 20 40 60 20 40 60 

x 2 + y 2 (m 2 ) 


FIG. 6: Vertex distribution of atmospheric neutrino candidates in the 
energy range of 30—1000 MeV for (a) pointlike events (\t < 1.2) 
and (b) short-tracklike events (1.2 < \t < 5.0). 


for M pulses) and A is the mean of the log(Ax 2 ) distri¬ 
butions evaluated with the p — >• VK + MC as a function of 
AT. logfXshape) is required to be more than —0.5. 

Criterion 2 is essential for removing the background from 
OD untagged cosmic-ray muons. These events concentrate in 
the top and equator regions of the detector because the water 
Cherenkov active volume outside these regions is relatively 
thin. Muons traversing these regions also tend to have short 
paths. Figure0shows the vertex distributions of atmospheric 
neutrino candidates in the energy range of 30—1000 MeV. The 
requirement of a pointlike hit timing distribution is effective at 
removing this OD untagged muon background within the full 
KamLAND volume. To maintain efficiency for K + decays at 
larger AT this fit must be relaxed, but this can only be done 
within the tighter fiducial volume defined by (2). 

Figure 0 shows the distributions of best fit energies of 
the two pulses (£ 1 , £ 2 ) for data, atmospheric neutrino MC 
events, and p -A VK + MC events, after all cuts except (4). 
The main background sources are misidentified single-pulse 
events from atmospheric neutrinos due to the limited multi¬ 
pulse separation power in the fit. The multipulse fit result 
for p —? VK + MC events (Fig. 0 right panel) should repro¬ 
duce the kinetic energy distribution for K + (£ 1 ) illustrated 
in Fig. 0 and for decay daughters of the K + ( £ 2 ). In the 
£2 distribution, there are two bands corresponding to stopped 
K + decays, K + —> p + v ll with a branching ratio of 63.55% 
(the kinetic energy of the p + is 152 MeV) and K + — > 7r + 7r° 
with a branching ratio of 20.66% (the sum of kinetic energies 
of the 7 r + and the 7r°-decay-induced 7 rays is 354 MeV). The 
difference between the visible energies and the real energies 
is explained by the nonlinear scintillation response, primarily 
Birk’s quenching. The anticorrelation between £1 and £2 is 
due to the multipulse separation uncertainty for the case of the 
smaller time difference between the first and second pulses or 
to contributions from other K + decay modes. Figure0shows 
the distribution of the shape parameter for data, atmospheric 
neutrino MC, and p —f VK + MC, after all cuts except (5). 
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FIG. 7: Distribution of best fit energies of the two pulses E\ and E 2 for data (left), atmospheric neutrino MC events (middle), and p — VK + 
MC events (right). E 1 provides an estimate of the K + kinetic energy, as illustrated in Fig. [I] and E 2 the kinetic energy sum for decay daughters 
of the K + , summarized in Tabled] The box shows the selection cuts for E\ and i? 2 - The detection efficiency and the expected number of 
background events were estimated from the MC event distributions for p —► VK + (right) and atmospheric neutrinos (middle), respectively. 


TABLE II: Detection efficiencies estimated from p —► uK + MC at 
each event reduction step for the five selection criteria. 



Period I 

Period II 

(1) Vertex cut ( x , y, z ) 

0.979 

0.978 

(2) Pointlike cut (xt) 

0.935 

0.932 

(3) Time cut (AT) 

0.658 

0.624 

(4) Energy cut (E u E 2 ) 

0.454 

0.447 

(5) Shape cut (T s h ape ) 

0.444 

0.432 


The atmospheric neutrino backgrounds from misidentification 
of single-pulse events tend to have lower L s h a pe as expected. 
The cut values for E\ and L s h a pe were chosen using MC only 
maximizing the sensitivity of the search with respect to the 
resulting variation of the detection efficiencies, the expected 
number of background events, and the expected statistical un¬ 
certainty. The expected number of atmospheric neutrino back¬ 
ground events is 0.9 ± 0.2 in total. No data remain after all 
cuts. 


V. EFFICIENCIES AND SYSTEMATIC UNCERTAINTIES 

Table |IT] shows the detection efficiencies estimated from 
p —> VK + MC at each event reduction step for the five se¬ 
lection criteria. The MC statistical uncertainty is less than 
0.5%. The largest inefficiencies, caused by the selection crite¬ 
ria (3) and (4), are highly dependent on the multipulse separa¬ 
tion power in the fit. The decrease of the detection efficiencies 
for the case of K + decay within a short time is illustrated in 
Fig. 0 The total efficiencies for p —> VK + in the LS are es¬ 
timated to be 0.444 ± 0.053 for period I and 0.432 ± 0.050 
for period II, considering the systematic uncertainties on the 
detection efficiency listed in Table Hill 

The number of protons is calculated based on the total LS 



lo §( L S ii ape ) 


FIG. 8: Shape parameter ( L s ha P e) distributions of data (black dot), 
atmospheric neutrino MC (thick red line), and p —> VK + MC (thin 
blue line) for (a) pointlike events (\t < 1-2) and (b) short-tracklike 
events (1.2 < Xt < 5.0). The contribution of the atmospheric 
neutrino background is normalized to the live time of data. 


volume (1171 ± 25 m 3 ) measured by flow meters during de¬ 
tector filling, the LS density (0.780 g/cm 3 ), and the hydrogen- 
to-carbon ratio (1.97) verified by elemental analysis. The 
nitrogen-to-carbon and oxygen-to-carbon ratios are negligi¬ 
bly small (< 6 x 10 -4 ). The 913 tons of LS mass contain 
7.74 x 10 31 free protons in hydrogen and 2.33 x 10 32 bound 
protons in carbon, and thus 3.11 x 10 32 protons in total. The 
uncertainty on the number of protons is 2.1%, determined 
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FIG. 9: Detection efficiency as a function of AT estimated from 
p —> vK + MC in period I. The decrease for AT < 7 ns, due to poor 
resolution of the multipulse fit, is the primary source of inefficiency 
for the p —> VK + search. 


from the LS volume measurement. 

The detection efficiency calculation from p —> VK + MC 
has associated systematic uncertainties. The MC code con¬ 
siders A' + -nucleon interactions in 12 C, which cause the K + 
loss with a probability of 2 . 1 % assuming the nuclear model 
given in Sec. lIII Al We conservatively include the size of this 
correction as a systematic uncertainty. The detector energy 
scale in the low-energy data is well calibrated by the radioac¬ 
tive source 7 rays, but in the energy scale of this analysis we 
must rely on the constraints provided by the fit to the spec¬ 
tral shape of atmospheric neutrino candidates up to 1 GeV 
shown in Fig. [3] In this fit we float the linear energy scale, the 
normalization factors, and the two systematic parameters pro¬ 
vided in the GENIE code: the axial mass for charged-current 
quasielastic scattering and the axial mass for neutral-current 
elastic scattering, which are inputs for the calculation of the 
nucleon form factor in nuclei. They are constrained to lie 
within their uncertainties. The energy scale factor is fitted to 
0 . 95^0 ^7 at ±lcr as indicated in Fig. [3] which corresponds to 
a 3.6%/1.9% uncertainty in the detection efficiency for Period 

m. 

The primary source of systematic uncertainty is a possible 
waveform deformation due to uncertainties in the scintillation 
process. The scintillation time profile is determined from the 
ionization density of charged particles interacting in the LS; 
thus the waveform shape depends on the type of particle in¬ 
teracting and its energy. As there are no waveform calibra¬ 
tion samples for K + scintillation, we computed the efficien¬ 
cies from p —> vK + MC for two extreme cases in which the 
scintillation time profile is constructed from the low-energy 
data for a 7 ray source and for a neutron source. The ion¬ 
ization density of a K + from proton decay is expected to be 
in between these two cases. We found a detection efficiency 
difference between these two of 11 . 1 % and assign this as a 
systematic uncertainty. The total systematic uncertainty on 
the lifetime calculation for p —*• VK + is 12.0%/11.6% for Pe¬ 
riod I/II, from adding the uncertainty sources summarized in 
Table Uni in quadrature. 


TABLE III: Estimated systematic uncertainties of the lifetime cal¬ 
culation for p —¥ VK + . The dominant uncertainty comes from the 
waveform modeling (11.1%) for the K + signal in the detection effi¬ 
ciency calculation. 


Source 

Period I 

Period II 

Number of protons 

2.1% 

2.1% 

Detection efficiency 

Nuclear model 

2.1% 

2.1% 

Energy scale 

3.6% 

1.9% 

Waveform model 

11.1% 

11.1% 

Total systematic uncertainty 

12.0% 

11.6% 


VI. RESULTS 

We find no event excess over background predictions in the 
p —» VK + search. The lower limit on the partial proton life¬ 
time is obtained from: 

t/B(p^VK + ) > T *i (3) 

where N p is the number of protons in the whole LS volume, 
Ti is the detector live time, e* is the detection efficiency for 
period I (i=l) and period II (i= 2) given in Table HT1 and numit 
is the estimated upper limit on the number of events for the 
p —> VK ' signal. Considering the large uncertainty of the 
background estimate, we conservatively calculate numit as¬ 
suming zero background. Using the Feldman-Cousins proce¬ 
dure |24l]. we calculate the confidence interval which accounts 
for the errors on nuisance parameters listed in TablelTIIl result¬ 
ing in numit = 2.44 events at 90% C.L. and the correspond¬ 
ing lifetime limit of r/ B(p —»• VK + ) > 5.4 x 10 32 years. 


VII. DISCUSSION 

We demonstrated a sensitive search for the proton decay 
mode p —>• VK + with a new search method using a large 
liquid scintillator detector. While the search sensitivity in 
KamLAND is about an order of magnitude smaller than a pre¬ 
vious search in Super-Kamiokande due to the more limited 
exposure, we confirm the conclusion that the SUSY SU(5) 
GUT with minimal assumptions predicting a partial lifetime 
of r/B(p —> vK + ) < 10 30 years is excluded. Extended 
SUSY-GUT models giving a partial lifetime of r/B(p —> 
77 K + ) > 10 33 years are still beyond our search range. How¬ 
ever, our search sensitivity has not been limited by back¬ 
grounds yet due to the effective background suppression of 
the coincidence method. It indicates that future increases in 
statistics are possible with LS detectors. 

Several future tens-of-kton LS detectors at different under¬ 
ground sites are proposed to investigate various physics top¬ 
ics. The LENA experiment was planned to perform a sensitive 
search for p —> VK + , and the potential of a LS detector with 
a fiducial mass of 50kton was investigated m. Assuming 
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a projected detection efficiency of 65%, the search sensitiv¬ 
ity is expected to reach r/B(p —> vK + ) > 4 x 10 34 years at 
90% C.L. within ten years [Tilt) , allowing the extended SUSY- 
GUT models to be tested. Similarly, the sensitivity for the 
planned JUNO experiment with a fiducial mass of 20kton 
LS is expected to reach t/B(j> —> VK + ) > 1.9 x 10 34 years 
at 90% C.L. within ten years ll25ll . Even if a de¬ 
tection efficiency of only 44% achieved in this study 
is applied to the searches in LENA and JUNO, the 
sensitivities are r/B(p —> VK + ) > 2.7 x 10 34 years and 
t/ B(p —> T7K + ) > 1.5 x 10 34 years at 90% C.L., respec¬ 
tively, which will still significantly impact models. The higher 
efficiency of 65% estimated in Ref. US is mainly due to 
the clear discrimination of atmospheric neutrino backgrounds 
based on pulse rise time measurements. From atmospheric 
neutrino MC studies, we find that the rise time stability in 
KamLAND is insufficient to provide efficient discrimination. 
In the cases of LENA and JUNO, kaon production induced by 
atmospheric neutrino interactions, discussed in Sec. II II Bl are 
non-negligible. However, if the background discrimination by 
counting the number of decay electrons suggested in Ref. Ill ill 
works, a background-free search will still be possible. 

VIII. CONCLUSION 

We have presented a search for the proton decay mode 
p —> vK + with the KamLAND liquid scintillator detec¬ 


tor. We found no evidence of proton decay and obtained a 
lower limit on the partial proton lifetime r/B(p —> VK + ) > 
5.4 x 10 32 years at 90% C.L. Due to KamLAND’s smaller 
detector size, the limit is less stringent than the latest Super- 
Kamiokande result, r/B(p —> vK + ) > 5.9 x 10 33 years 
at 90% C.L. However, the detection efficiency of 44% in 
KamLAND is about two times better than that of Super- 
Kamiokande for background-free searches, and the systematic 
uncertainty of 12% is almost free from nuclear modeling un¬ 
certainties. Based on these results, we expect future searches 
with tens-of-kton liquid scintillator detectors to be compet¬ 
itive with the water Cherenkov detectors. Such improved 
searches will test SUSY-GUT models more stringently. 
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